Fermi liquid identities for the Infinite U multichannel Anderson Model 
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We show how the electron gas methods of Luttinger, Ward and Nozieres can be applied to an 
SU(N) , multi-channel generalization of the infinite U Anderson impurity model within a Schwinger 
boson treatment. Working to all orders in a 1/N expansion, we show how the Friedel Langreth 
relationship, the Yamada-Yosida-Yoshimori and the Shiba-Korringa relations can be derived, under 
the assumption that the spinon and holon fields are gapped. One of the remarkable features of this 
treatment, is that the Landau amplitudes depend on the exchange of low energy virtual spinons and 
holons. We end the paper with a discussion on the extension of our approach to the lattice, where 
the spinon-holon gap is expected to close at a quantum critical point. 
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I. INTRODUCTION 

The standard model of interacting electronic systems 
is based on the theoretical framework established by 
Landau^ in the 1950's, whereby the properties of inter- 
acting fluids of fermions are linked to non-interacting 
fermions by adiabatically turning on the interactions. 
This is the basis of the Landau Fermi liquid, and it leads 
to the concept of quasiparticles. It also provides the for- 
mal basis for the diagrammatic description of interacting 
Fermi fluids. In the early 1960's, Luttinger and Ward 2 
showed how many aspects of the Landau Fermi liquid 
theory - such as the relationship between fermion density 
and Fermi surface volume, the relationship between the 
linear specific heat capacity and the quasiparticle density 
of states, could be deduced from first principle resum- 
mations of diagrammatic perturbation theory to infinite 
order. 

Many new electronic materials discovered over the past 
two decades, such as the high temperature superconduc- 
tors, colossal magneto-resistance and heavy electron ma- 
terials, do not fit so naturally into the Landau scheme. 
In these materials, the interactions between the fermions 
become so large they project out large tracts of the many 
body Hilbert space. For example, in many narrow band 
transition metal or f-electron materials, the charge fluc- 
tuations of the localized orbitals are restricted to certain 
valence states, and in local moment systems, the elec- 
trons in a localized moment are restricted to states of 
fixed definite occupancy. In these systems, the valid- 
ity of the adiabatic approach is increasingly in question, 
and a more appropriate starting point is one where the 
Coulomb interaction has been taken to infinity from the 
outset, explicitly removing certain states from the Hilbert 
space. For example, in the infinite U Anderson model, 
states of double occupancy are excluded whereas in the 
Kondo model, all charge fluctuations of the local moment 
are excluded, as illustrated in Fig. 1. 

When cast as a field theory, infinite U models be- 
come gauge theories in which the Hubbard operators 
that create and destroy the highly correlated electrons 



FIG. 1: Schematic diagram illustrating the evolution of the 
impurity f-spectral function with the strength of the repul- 
sive interaction U. Shaded ares denote the regions where 
the impurity f-spectral function is large. As U is increased, 
high-energy regions of the Hilbert space are projected out of 
the low-energy model, giving rise to the infinite-!/ Anderson 
model, and Kondo model as effective low energy theories. So 
long as adiabaticity is preserved, the spectral weight at the 
Fermi surface is preserved, no matter how large U becomes. 



are written as composite products of "slave operators" . 
In such an approach, the constraints on charge fluctua- 
tions become conserved local quantities associated with 
local gauge invariances and the basic fields entering into 
the Hamiltonian are "fractionalized" fields, which carry 
either spin or charge, but not both. The gauge theory 
approach to strongly correlated electron systems poses a 
number of important technical and qualitative challenges. 
In particular, how does the Landau Fermi liquid emerge 
within a Gauge theory of unbroken symmetry, and what 
happens to the spectrum of gauge particle excitations ? 
This question assumes an increased importance in the 
context of recent speculations that fractional particles 
may become free at a quantum critical point£ii£i£. 

The last ten years have seen considerable progress 
in methods that combine the Luttinger Ward approach 
with a new generation of large N expansions which use 
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Schwinger bosons, rather than Abrikosov pseudo-fermions 
to represent the localized moments^^ ' 10 ! 11 . We now 
briefly review some of these recent developments. Cox 
and Ruckenstein were the first to apply a large TV ap- 
proxach to a multi-channel Kondo single impurity Kondo 
model, usinga fermionic spin represenation. Parcollct 
and Georges Q developed the original formalism for the 
Schwinger Boson large N approach to this model, em- 
ploying a multi-channel Kondo model, focusing their in- 
terest mainly on the overscreened case @ , built upon this 
earlier work to develop a Luttinger Ward approach, show- 
ing how it can be used to derive a Friedel sum rule for 
the single impurity model and a Luttinger sum rule for 
the lattice, assuming that the Fermi liquid state could be 
treated in a 1/N expansion. Rech et al. [l(| showed how 
the fully screened Kondo impurity model could indeed be 
treated in the large N limit using the Schwinger boson 
method, demonstrating the emergence of Fermi liquid be- 
havior. Lenanon et al. (Tlj extended these methods to 
the infinite U Anderson model, showing how a Baym 
Kadanoff approximation scheme can be developed from 
the 1/N expansion of the Luttinger Ward functional, us- 
ing the leading order term to develop a conserving ap- 
proximation at finite N. 



One of the key insights emerging from the large N 
solutions to the fully screened Kondo and Anderson im- 
purity models, is that the development of a Fermi liquid 
is accompanied by the formation of a gap A g of order the 
Kondo temperature Tr-, in the spectrum of the fractional 
particle s 10 ! 11 . 



In this paper we explore the general consequences of 
this observation for the infinite U Anderson impurity 
model. In particular, we show that the the assumption 
of a gap in the fractionalized particle spectrum enables 
us to extend the Luttinger Ward approach beyond the 
leading orders considered so far, to all orders in the 1/N 
expansion. In our approach, the expansion in powers of 
1 /N plays a completely analogous role to the expansion 
in powers of the interaction. The Luttinger Ward ap- 
proach was first applied to the finite U Anderson model 
by Yamada, Yosida, Shiba and Yoshimori ) 12 i 13 i 14 who 
used conservation laws to develop a set of conserving 
identities that apply to all orders in the strength of the 
interaction U. Here, taking advantage of these new in- 
sights, we show how a parallel set of results can be ob- 
tained for the infinite U Anderson model. 



Many structural aspects of our work, including the 
interplay between virtual fractional excitations and the 
composite heavy quasiparticles may enjoy a wider appli- 
cation to the Anderson and Kondo lattices. These are 
points that we return to in a discussion at the end of the 
paper. 



II. SUMMARY OF KEY ASPECTS OF THE 
PAPER. 

The starting point for our work is the infinite U Ander- 
son model, which we formulate by representing the Hub- 
bard operators as a product of slave fermion ("holon") 
and a Schwinger boson ("spinon" ) 15 ' 16 X a o = b\ t \ where 

creates a "spinon" with spin component a and x cre- 
ates a charged "holon" . Since there is no natural per- 
turbative expansion in terms of the parameters of the 
original model, we consider a family infinite U impurity 
Anderson models with SU (N) spin symmetry, employing 
a 1/N expansion around the large N limit. There are a 
number of special tricks that need to be carried out in 
order to realize such a large N expansion 

• The spin S of the moment, described by the number 
of Schwinger bosons rib = 25 must be allowed to 
grow with N, so that 2S/N remains finite. 

• More surprisingly, the number of electron scatter- 
ing channels K must also grow with A 8 , in order 
to screen the large moment 5. The special case 
25 = K describes the fully screened model 1 ^. 

The introduction of K channels requires that we gener- 
alize the Hubbard operators as follows 1 ^ 



Xv 



(1) 



where v G [LA] is the channel index. X av creates a 
localized electron with spin component a and channel 
index v. 

We formulate these features into an infinite U Ander- 
son impurity modelii as follows, 



H = 



(2) 



Here ct creates a conduction electron with momentum 

k, channel index v G [1,A], and spin index a € [l,iV]. 
The term 



H r 



kva. 



H.c 



describes the interaction of the conduction electrons with 
a magnetic ion located at the origin. The y/N de- 
nominator in the hybridization ensures a well-defined 
large N limit. The energy of a singly occupied impu- 
rity eo is taken to be negative. The conserved operator 
Q = J2 a blfia + xtXv generalizes the no-double oc- 
cupancy constraint of the infinite U Anderson model, by 
restricting the "valence" nj, < Q. Q is also the maximum 
size 5 = Q/2 of the local moment that can develop at 
each site. The condition Q — K is required for perfect 
screening of the local moment. 

The spinon and holon operators carry the conserved 
charge Q, and in order to discuss these fields, we need to 
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consider the enlarged Fock space of general Q ^ K. The 
partition function for the general problem is given by 



J Q 



Tr [V Q e 



-PHI 



where the projection operator Vq imposes the constraint 
Q = Q. In so doing, we are really considering a whole 
family of Anderson/Kondo models which includes under- 
screened models where Q > K is larger than the number 
of screening channels and the overscreened models where 
the number of channels K > Q exceeds the number of 
spinons per site. 

We shall make the key assumption that the fully 
screened state where Q — K develops an additional 
stability with respect to the over and under screened 
states where Q =/= K. This assumption is motivated 
by the discovery of a spinon/holon gap in the large N 
Schwinger boson solution of the one and two impurity 
Kondo models^ and the leading Baym-Kadanoff approx- 
imation to the infinite U Anderson model 11 . The struc- 
ture of these solutions lead us to believe that the gap in 
the fractional particle excitation spectrum is robust in its 
extension to finite N. The emergence of a gap in the infi- 
nite N limit provides a vital infra-red cuttoff to the fluc- 
tuations of the fractional excitations. This cut off plays 
a dual role: it stabilizes a Fermi liquid with a large Fermi 
surface, whilst at the same time, providing a vital low- 
energy cut-off to the fluctuations of the fractional holon 
and spinon fields removing them from the low-energy 
Hilbert space. The presence of the gap guarantees that 
the 1/N expansion contains no non-perturbative infra- 
red divergences, and this in turn guarantees its survival 
over a finite range of 1/N values. 

Suppose E is the energy of the ground-state with 
Q = K and E± be the ground-state energies of the state 
where where Q — K ± 1 is modified by one unit, then the 
"commensuration gap" associated with complete screen- 
ing of the moment at site jo is 



1 



A g = -[E + + E_ - 2E ]. 



(3) 



This gap is manifested as a gap in the spectral functions 
of the spinon and the holon as illustrated in Fig [2j With 
a gap we can calculate the low-temperature properties of 
the model in the grand-canonical ensemble where remove 
the constraint and associate a chemical potential A with 
the conserved charge Q 



Z a [X] = e-^W = Tr 



(4) 



Although the chemical potential A only imposes Q — 
Q on the average, the existence of commensuration gap 
guarantees that there is a finite range of A 

A e [E — E-,E+ — E ], 

for which the constraint is imposed with exponential ac- 
curacy once T« A. 
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FIG. 2: Illustrating the spectral gap in the holon and spinon 
spectra. Creation of a holon or spinon increases Q — » K + 1, 
whereas the destruction of a holon or spinon decreases Q — > 
K — 1. 



Our approach makes use of Coleman-Paul-Rech gen- 
eralization of the Luttinger Ward Free energy func- 
tional, to the gauge theory of interacting bosons and 
fermion o 10 i 18 ' 19 which can be compactly written as 



F = TStr [In (-G' 1 ) + E0] + Y[0\, 



(5) 



where Str [A] = Tr[As] — TtL4f] denotes the supertrace 
of a matrix containing both bosonic (B) and fermionic 
(F) components (where the underline notation is used 
to denote a sum over internal frequencies and a trace 
over the internal quantum numbers of the matrix) . Q = 
{G^ 1 — S) _1 is the matrix describing the fully dressed 
Green's function of all elementary particles and fields 
entering the Lagrangian, including the slave particles, 
where E is the self-energy matrix and Go the bare prop- 
agator of the fields. Diagrammatically the propagators 
for the Anderson model are denoted as follows 



G r 



lUJr, 



G b (ih> n ) 

GJiu n ) 



■->-- 



In the above expressions we display two alternative nota- 
tions for the propagators. In the first column, we adopt a 
traditional notation, using continuous, wavy and dashed 
lines to represent the conduction, spinon and holon prop- 
agators respectively. In the second column, we display 
the alternative "railway track" notation that we use in 
cases where it is necessary to clearly show the flow of 
spin and charge in the respective propagators. At the 
interaction vertices 



V 



N 



V 



N 



the spin and charge of the electron fields divides up 
amongst the holon and spinon fields. In the Feynman di- 
agrams, each vertex is associated with a factor iV/y/N. 
The generating functional Y[Q] may be written as the 
sum of all closed-loop two-particle irreducible skeleton 
Feynman diagrams. These can be ordered in a 1/N ex- 
pansion as follows 




Each blue loop over spin gives a factor of N, whereas 
each red loop over charge gives a factor of K = kN. 

The functional Y[Q] is the generating functional for 
the self-energies of the fields, so that variation of Y[Q] 
with respect to the full Green's functions Qq of the con- 
duction, spinon and holon fields, = —r]^/3 SY/6G(, 
self-consistently determines the self-energies of these 
fields&i^, where rj^ = ±1 for bosons/fermions. We refer 
to Y[Q] as the Luttinger Ward functional. Approxima- 
tions to this functional are the basis of Kadanoff-Baym 
conserving many body approaches 2 ^. 

The key difference between our new approach, and 
the classic Luttinger Ward approach to the interacting 
electron gas, is the appearance of the holon and spinon 
fields. In the Luttinger Ward approach, the Y[Q] is a 
generating functional for the self-energies of the fields. 
The second-derivative of the generating functional 
with respect to the Green-functions generates the 
off-shell scattering vertex function between the fields 
A = f3 2 5 2 Y[Q}/SQ 2 . Nozieres and Luttinger— demon- 
strated that the multiple forward scattering generated 
by this vertex gives rise to the on-shell Fermi liquid 
amplitudes between the electron quasiparticles. The 
Landau amplitudes are determined as solutions of the 
Dyson equation represented by the Feynman diagrams, 

® -® + £03 

This feature is preserved in the new approach, but the 
off-shell scattering vertex also involves virtual holon or 
spinon pairs, as shown below: 
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In this way, scattering off the virtual fractionalized exci- 
tations determines the scattering amplitudes of low en- 
ergy Fermi liquid. Remarkably, despite these effects, all 
of the key conservation properties of the Landau Fermi 
liquid are preserved, so long as the fractionalized excita- 
tions remain gapped. 

In the context of impurity models, the Landau Fermi 
liquid becomes Nozieres' "local" Fermi liqui d 22 i 23 , in 
which the energies of each electronic quasiparticles are 
expressed in terms of the elastic scattering phase shift 
5 ua , where v is the channel index and a the spin index. 
If Ae is the spacing of states in the continuum, then the 
energy of each electron in the continuum is shifted down 

by an amount — ^ s "™' i J Ae. The Fermi liquid theory is 

then defined in terms of the dependence of the phase 
shifts on the quasiparticle occupancies 

Sva[e, {nkau}] = + e6' + ^ <j>va,v'a'&rikui<x> , (7) 

where the impurity Landau parameter 

<t>ua,u>u' = (fj a; v'a'\cj>\i>, a; v'a) (8) 

is the expectation value of the interaction energy between 
quasiparticles of flavor/spin (y,a) and (z/,a'). 

By comparing the Free energy of the infinite U Ander- 
son model obtained in the Luttinger Ward approach, we 
are able to diagrammatically identify each of the terms 
in the above expansion. In particular, the conduction 
electron phase shift is given by 

S va (e) = Imln(l - i7rpE c (e - i0 + )), 

where E c (e) is the conduction electron self-energy, p the 
density of conduction electron states at the Fermi sur- 
face and + is a small positive infinitesimal. Wc are also 
able to identify the Landau parameter with the on-shell 
forward scattering amplitude between quasiparticles 

<Pua,u'a' — F va y a >(u} = 0, U) = 0). 

As in the finite U Anderson model, each of these quan- 
tities enjoys expression entirely in terms of conduction 
electron states. 

The ability to link the quasiparticle physics with the 
gauge theory diagrammatics enables us to derive the key 
conservation relationships for the infinite U Anderson im- 
purity model, working to all orders in the 1/N expan- 
sion. These relationships survive in the leading Kadanoff 
Baym 11 approximation that defines the large N limit, 
and we can illustrate them explicitly in practical calcu- 
lations. The key conservation laws that appear from this 
approach are 
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The Friedel sum rule^i^^ 

7r K — n 



N K 



\ND 



(9) 



III. PRELIMINARIES. 



Our starting point is the infinite U Anderson model. 



• The Langreth relation^ 

Imt(e - i0+)| e=0 = sin 2 <5 c /7rp, (10) 
where p is the density of states of the electron fluid. 

• The Yamada-Yosida-Yoshimori relationship 12 ! 13 
between the spin, charge and channel susceptibili- 
ties and the linear specific heat coefficient 



n 



E6-C-. C — 
k kvoi kvot 

kva. 



V 



N ^ 

kvGL 



[4 x% 



H.c 



bib a -X{Q-2S), 



(12) 



NK 



7 



A' 



N- 



K 



1 Xf . Xc 



where the spin components a are taken to run from —j 
to +j, where N = 2j + 1, and the flavor components v 
are taken to run from — / to +/, where K = 2f + 1. 
The bare conduction electron, spinon and holon Green 
functions are given by 



7° N + K X °s K + N X ° f X° c 



(11) 



*2*W = vh: 

,(0) 



z - (eo - A) 
1 



z + X 



(13) 



where we use the symbol z to denote the analytic exten- 
sion of the Matsubara frequencies into the complex plane. 
In an impurity model, it is generally convenient to trace 
the conduction electron propagator over its momentum 
indices to obtain the local propagator. The bare local 
conduction propagator is given by 



i 



• The Korringa Shiba relationship^ between the dy- 
namical spin susceptibility and the impurity sus- 
ceptibility x"M/wL =0 = (Nn/K)( X /N) 2 . 

The structure of the paper is as follows. In section (3) 
we set up various formal preliminaries. In (4) we reca- 
pitulate the relationship between conservation laws and 
Ward identities, using them to rederive a general set of 
Friedel sum rules. In (5) we show how these can be used 
to derive the Langreth 2 ^ relation for the scattering t- 
matrix and a general set of Friedel sum rules. In (6), we 
introduce the off-shell interaction vertices amongst the 
electrons and the fractional particles. Using the Friedel 
relations we express the charge, spin and flavor suscepti- 
bilities in terms of these vertices. In (7) we use our gap 
hypothesis to derive an expression for the low tempera- 
ture Free energy and use this to relate the specific heat, 
with the spin, charge and flavor susceptibilities. In (8) we 
identify the Nozieres interaction parameters 2 ^ 3 , with the 
on-shell interaction vertices of the gapless electrons and 
use our results to derive an identity between the linear 
specific heat coefficient and a linear combination of the 
spin, charge and flavor susceptibilities. Our result here 
is a generalization of the Yamada-Yosida-Yoshimor i 12 i 13 
identities for the finite U Anderson model. In section 
(9) we derive the Korringa-Shiba relationship 1 '! between 
the uniform susceptibilities and the power-spectrum at 
low energies. Finally in (10) we discuss the extension of 
these ideas to a lattice environment. 

I 

Each of these self-energies is obtained by differentiating the Luttinger Ward functional, shown in ([6]) with respect 
to the Green's functions 



=Fz7rp forz = w±iO + , (14) 



where the expression on the right is the large band width 
limit of the local propagator, choosing the negative sign 
for the retarded propagator (z above the real axis) and 
the positive sign for the advanced propagator [z below 
the real axis). 

When interactions are turned on, each of the fields 
acquires a self-energy correction in the Full propagators 



g c (z) = [[gpHz)]- 1 -^)]- 1 , 

g b (z) = [z-(6 -A)-£ fc (z)]-\ 

g x (z) = [z + \-x x (z)}-\ 



(15) 



E&M = -0 



E x (u) = +[3 



SY 



sgb(u) 



SY 

sgju) 





0(1) 



0(K/N 2 ) 



where the cross indicates the line which is eliminated by the functional differential. Although our discussion today 
will focus on keeping all orders in this expansion, it is interesting to briefly reflect on the leading terms in the 
diagrammatic expansion. The leading order contributions to the self energies derive from the virtual separation of the 
electron into spinon-holon pairs. These terms are already sufficient to develop a conserving approximation with all 
the key features of the infinite order resummation. The second-order terms in this expansion contain the sub-leading 
effects of electron-electron scattering. 



IV. FRIEDEL SUM RULES AND LUTTINGER 
WARD IDENTITIES. 

The Fermi liquid description of this model centers 
around the relationship between the conduction electron 
phase shift and the thermodynamics, and the charge, spin 
and flavor susceptibilities. We begin with a derivation of 
the Friedcl sum rule. This part of the paper closely fol- 
lows earlier work on the Kondo model 

In addition to the local conserved charge Q = rib + n x , 
there are three globally conserved, physical quantities: 
charge (Q c ), magnetization (M = Q s ) and flavor (F = 
Q F ) associated with this model, which we write as 



kva v 
Q S = E 5 C L a C k,a + E &b ° b <* = 

Q F = E* c L c ^-E^x„ (it) 



kva 



respectively, where a = sign(a) is one for "up spins" 
and minus one for "down spins" , and v = sign(f) is one 
for flavors v > and minus one for flavor index v < 0. 
Each of these conserved quantities is associated with a 
corresponding external field, which we introduce into the 
Hamiltonian by writing 



H = U- A f iQ c - BM - B F Q F , 
= H- £ B A Q A , 

A=C,S,F 



(18) 



where we introduce the notation (B C \B S ,B F ) = 
(A/x, B,B F ). We shall write each of these conserved 
quantities using the short-hand 



(19) 



where ^ = (ct, tf,x^) denotes the complete spinor of 
electron, spinon and holon fields. We can write the ex- 
pectation value of these operators in terms of the triplet 





conserved quantity 


particle 


q c 


<7 S 


F 

1 


c 


1 


Q 


V 


b 





a 





X 


-1 





— V 



TABLE I: The diagonal components of the conserved charge, 
spin and flavor operators q c , q' s and q F ', resolved into their 
conduction (c), spinon (b) and holon (x) components. 



of electron, spinon and holon Green functions 

g = {g cU „g b ,g x }, m 

as 

(Q A ) - Tr {^q A i>g} = -TStr [q A g] , (21) 

where g is the thermal density matrix and the supertrace 
denotes a trace over each particle species (with a minus 
sign for fermions) and their quantum numbers, and the 
underline beneath the super trace denotes a summation 
over Matsubara frequencies, 

= ^ Tr [q A g B ^n)] - ^ Tr [q A g F {iuJ n )\ . 

An alternative, and convenient way to formulate the 
conserved quantities at absolute zero, is to replace q A by 
its more general form, q A — > q A j(uj) where 



so that 



{Q A )=-TStr[q A 1 g]. 



(22) 



(23) 



This is a particularly useful device in the single impurity 
model where we work with the local conduction electron 
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propagator. In this case the charges do not couple triv- 
ially in the local propagator and we need to introduce 
frequency dependent vertices with the external fields. To 
see this, note that the bare local propagator becomes 



g° c ( z ,B) = £[( 



z - e k 



AdAi-1 



= g° c (z) - [g° c (z)f(B A q A 1 {z)) + 0(B 2 ), 



where 



lc{z) 



dz 



(24) 



is identified as the frequency dependent vertex function. 
The conduction electron vertex 7 c (w) vanishes in the 
wide band limit j c — > 0. This is the basis of the fa- 
mous "Anderson Clogston" compensation theore m 27 ' 28 , 
according to which the conduction band polarization of 
charge, spin or flavor degrees of freedom in the ground- 
state 



(Qc 



Tr Iq^Gc 



0(T K /D), 



where D is the electron band- width and Tj( is the Kondo 
temperature. 

Ward identities are properties of the zero temperature 
self-energies, Green's functions and scattering vertices 
that result from the the conservation laws. A key step in 
making the transition from finite, to zero temperature, is 
the replacement of Matsubara sums at finite temperature 
by continuous integrals along the imaginary axis at low 
temperature 



A(ia n 



da 

du> 
2ttI 



A{ia) 



(25) 



Here a n represents the Bose Matsubara frequencies v n — 
2nirkBT along bosonic propagator lines, or the Fermion 
Matsubara frequencies w„ = (2n + l)7rfcsT along fermion 
lines. Our ability to make this replacement presumes the 
existence of well-defined scales associated with the ex- 
citation spectrum. This is a situation that may not be 
satisfied at a quantum critical point, or an overscreened 
Kondo model where the spinon and holon propagators 
may exhibit E/T scaling, with temperature as the rel- 
evant excitation scale. However, away from a quantum 
critical point, the conduction electron lines are character- 
ized by a well-defined Fermi energy and the existence of a 
well-defined gap in the spinon and holon spectrum guar- 
antees that the continuous version of the Matsubara sums 
is valid for these lines too. Provided the replacement l[2"5"]) 
is allowed the zero temperature conserved quantities may 
be written in the form 



(Q A ) = 



/ioo j 



(26) 



The starting point for the development of Ward iden- 
tities is the Luttinger Ward functional. The conserved 
charges are conserved at each vertex in the diagrams con- 
tributing to the generating function Y. For the quantity 
Q A , the charge q A is associated with each propagator line 
£7f inside the functional. When we shift the frequency of 
each propagator by a constant differential times the cor- 
responding charge 



Slu a = Suj qf, 



(27) 



the generating functional is unchanged. Since the deriva- 
tive of the generating functional with respect to the 
Green's functions is the self-energy, this leads to the 
Ward identity 



SY[G] 
Slu a 



and integration by parts then yields 



(28) 



(29) 



Equations (|26p and ([29)) sum up to an integral over a full 
differential, rotating the integration contour around the 
negative real axis we get the Friedel sum rules 



400 duj „. f, A d£ 



1 



(Q A ) — Im Str{g A ln [-g-\-M+)] } 



(30) 



In impurity problems, we are interested in change in these 
quantities that results from coupling the impurity to its 
environment, which is given by 



(Q A ) 



= — Im Strain [S-7So _1 ]}U- 



i0+ 



= -- Im Str{g A ln[l - g (uj)T,(uj)]} 



(31) 



The quantities on the right hand side are the "phase 
shifts" of the electrons, spinons and holons. If we identify 



5 c = Im{m[l-Soe(a;)£fM]}L = _ i0+ 
then we can write 



(C = c,b, X ) 
(32) 



t — ' 7T — ' 7T » — ' 7T 

At first sight, this is very different from the expression 
we would expect for a Fermi liquid, where there are no 
holon or spinon contributions. However, so long as the 
spinons and holons are gapped, their corresponding phase 
shifts arc cither n or and do not vary with the external 
fields. General arguments lead us to believe that in the 
impurity model, the spinon phase shifts are zero, while 
the x phase shifts are tt. To see this, consider the Friedel 
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sum rule for the conserved charge Q 
{ie %, q x ) = (0> 1) !)■ In this case, 



rib + n x , for which 



AQ = 2S = K 



7T 7T 



(33) 



Provided all the phase shifts are positive definite, then 
this quantity can never exceed K and the maximum value 
is only attained if 6b = and 5 X — n. If the spinon phase 
shifts are zero, then for a perfectly screened impurity, 
with K = 2S it follows that the holon phase shifts equal 
7r and the spinon phase shifts vanish so long as the spinon- 
holon gap is preserved. 

The general sum rules then become 

AO c = \ ^ cua - \ = \ ^ cva — K 



q A 



^CVOL 



(A = S,F), 



where the conduction electron phase shift is 
5 cva = Imln[l - G°T, cl/a ] 



\u=—%0+ ■ 



(34) 



(35) 



The x fermion contribution only enters into the charge 
Friedel sum rule, since Yl v = an d J2 V Qx'i^ = ^ - 
The first sum rule expresses the conventional Friedel sum 
rule. In general, the impurity charge AQ C contains a 
conduction electron and a holon contribution 



AQ = AQ c e 



(36) 



According to the Anderson Clogston theore m 27 ' 28 , the 
first term is of order Tk/D, where D is the band- width, 
and vanishes in the infinite band- width limit. The re- 
maining term is simply the charge associated with the 



empty /— states at the impurity, AQ = — 
ing the phase shift expression for AQ C in eq. 
Friedel sum rule can then be rewritten as 

Eir =*-<»*>• 



Us- 
the 



(37) 



In the case where there is no magnetic or flavor po- 
larization, all phase shifts are equal, 8 cva = S c and 
Si/a ficva = NKS C , so that the Friedel sum rule becomes 



6 C = - 
N 



1 - 



K 



(38) 



V. LANGRETH SUM RULE AND THE KONDO 
RESONANCE 

The emergence of the Kondo resonance is one of the 
signature features of the Kondo effect. The "Langreth 
sum rule"- 6 , links the Friedel sum rule to the formation of 
the Kondo resonance by providing a rigorous constraint 
on the size of the spectral function associated with the 
localized state. 



The t- matrix of the conduction electrons is determined 
by the spectral function of the electrons localized in the 
magnetic moment, 



' N 







(TXi„(T)X a „(0))e^ T dT (39) 



We shall assume here that there is no magnetic or flavor 
polarization, so that the t-matrix is identical for all v and 
a. The t-matrix is determined diagrammatically, by the 
process of repeated scattering off the localized state, and 
is given by 



(40) 



or 



t(z) = E c (z) + E c (z)G(°)(z)S c (z) + 
= S c (z) 
1 - G°(z)E c (z) 



(41) 



From p5[) we see that denominator in the t-matrix is also 
the argument of the scattering phase shift. In the limit of 
a broad band- width, we can replace g£ ^ (u> + iS) = —iitp, 
so that the phase shift and t-matrix become 



t(uj + ie) — 



T, c (oj + ie) 



1 + mpY, c {ij + ie) 

Im ln[l + i7rpE c ] | aJ= _ i0 + • 



(42) 



If we examine the decay processes associated with the 
conduction electron self-energy, we see there are two 
types of process. In one type, the electron decays into 
holon, spinon combinations, but since these are gapped 
excitations, they produce no contribution to the inelas- 
tic electron decay processes at the Fermi energy. The 
only remaining decay processes involve the production of 
electron- hole pairs, such as 




Ml -/*)(!-/*") 
(43) 



where the internal loop of gapped fractional quasipar- 
ticles does not contribute to the low-energy decay pro- 
cesses. These are the well-known decay processes of the 
conventional Fermi liquid, producing an inelastic decay 
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rate that grows quadratically with energy and temper- 
ature £"(w) ~ I^ 2 + ( 7r '^) 2 ]- At zero temperature and 
frequency, the conduction electron self-energy is thus en- 
tirely real, S' c '(0) = 0. So from gl]), tan<5 c = irp£ c (0), 
and the on-shell t-matrix is given by 



t(0 + ie) = 



tan<5 c 



sin S c 



■np (1 + itan(5 c ) irp 



and the imaginary part of the t-matrix is then 

sin 2 S r 



lmi(0 + ie) = 



TTf) 



(44) 



(45) 



the above expressions with respect to the corresponding 
fields X A = d(Q A ) imp /dB A at zero field B A = 0. Since 
the spinon and holon are gapped, their phase shift is 
pinned to the constant values of and 7r, and therefore 
the only contribution to the susceptibilities comes from 
derivatives of the conduction electron phase shift S = 
Imln[l-G°£ c ] : 



XA 



-Yq A 

TT ' * 



dS al , 



dB A 



(46) 



B A =0 



The existence of this quantity as an adiabatic invariant 
is well-known in the finite U Anderson model, and its 
appearance in the limit that U — ► oo is an important 
test of these methods. 



VI. STATIC SUSCEPTIBILITIES 

The susceptibilities of the impurity contribution to the 
conserved charges can now be obtained by differentiating 

I 



Explicitly the derivatives of the phase shift are given by 

, (47) 



dS =- l ^ c §^-im(g c + g c ) lc q A 



dB A 



-i0+ 



and in the wide bandwidth limit in where 7 C vanishes 



dS 
dB^ 



-ImGc 



d^r 



dB A 



(48) 



-i0+ 



To calculate the derivatives of the self energies it is useful to introduce the bare and full off-shell interaction vertices; 
The bare interaction vertex A is defined as the functional Hesseyan of the LW functional 



A aa /(cj,a/) = (2Tri) 2 r] a ri a , 



6 2 Y[G] 



SGa'(to') 



(49) 



where a and a' list the particle species c, b and x an( i their corresponding indices a and/or v, while (r]b,r] c ,r) x ) — 
(1,-1,-1) grade the anticommuting character of the Fermi fields. The 2ni factors come to compensate for the 
integration measure and replace the 1/T factors in the zero temperature limit when the discrete Matsubara frequencies 
are replaced by the continuous imaginary variable lo. For example the interaction vertex between an electron and a 
spinon is 



A c6 (w, J) = 



5g c {cj)Sg b (cu') 



0(1/N) 



0(K/N 3 ) 



0(KSN 3 ) 



0(K/N 4 ) 



0(K/N 4 ) 



The derivatives of the self energies with respect to the fields B A are expressed in terms of the bare interaction 
vertex by 



d£g(u>) 
dB A 



dB 1 



(50) 



adopting a heuristic matrix notation we may write this or 



dB A ' 

which implies that 



(51) 



(52) 



where 



~d~B A 



= -Trjg 2 ( iq A ) 



r= (i - A v g 2 ) 'a 



(53) 



(54) 
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is the solution to the Dyson equation 

T = A + Ar]g 2 T. 
Restoring the indices, we then have 

T aa ,(w,w') J ^-A QQ »(w,w") x 

a" 

{2m «J „ i0 , 8(lj" - oj') + r,a»Ql»{u") r o „ ,(w",w')} . 

It is instructive to represent this equation diagrammati- 
cally. For example, the vertex with incoming and outgo- 
ing electrons is r cc , represented as 

aval a'v'ai' ava> a'V'ef ava> a"m" ctVctf 




v'to" a'v'oj" avw 



a"v"co" a'v'ay 



The self energy derivatives are then reduced to 
dZ a (uj) ^ f du' , 2 , , , 

a' 

(55) 

Unlike the Fermi liquid in the original works of Luttinger 
and Ward, the on-shcll interaction vertex which repre- 
sents the Fermi liquid is normalized by the appearance 
of intermediate spinon and holon states. 

Incorporating the expression for the self energy deriva- 
tives within Eqs. (|48|) and (|46j). we can write the suscep- 
tibilities as 

A ImC c f diJ ^ — A A \ s-r2 / \ A 

olvci 

(56) 

Since we work in the wide band limit and 7 C vanishes 
the summation over 7 a q a gives non-vanishing contribu- 
tion only for one particle species 

o ImG, 



-i^Vfis' / ^r CQ . fcQ ,(o,^ b 2 H,(57) 

7T J 27TZ 

cxa' 



ImG, 



ImG, 



^ / 



doj 
2ni 



T cv , xv ,{Q,u)gl(u). (59) 



temperature Green's function Q = 0\t=o > the lead- 
ing change in the Free energy at finite temperatures 
AF = F(T) - F(T = 0) can be computed solely from 
the finite temperature trace of the Logarithm of the zero 
temperature Green's function 



AF = 



TStr {in 



T 

T=0 



(60) 



This is a quite non-trivial result, and for clarity, here 
we repeat the original Luttinger Ward derivation in the 
context of the current work. At small finite temperature 
we have two things to keep track of in the free energy, the 
discreteness of the frequency summations, and the tem- 
perature dependence of the Green functions. However, 
since the Free energy is stationary at low temperatures 
with respect to small changes in Q, to calculate the low 
temperature corrections to the Free energy, we may use 
the zero temperature Green function Q = Q\t=o in the 
Luttinger Ward expression 



F[T~0] = rStr hx[-g 



where (Q 1 



tg +Y[g,T], (6i) 



* — g ) = S is the zero temperature self- 
energy. 

Small finite temperatures introduce a discreteness into 
the frequency summations that produces corrections of 
order T 2 in the Free energy. To keep track of these 
changes we must write 

/ioo j 
— [1 + eHJAM, 



where 



e(ui) = {2niT^25{iuj + a n ) - 1) 



(62) 



keeps track of the discreteness. Luttinger and Ward made 
the insightful observation that at low temperatures, one 
can regard e(ui) to be an infinitesimal, so that inside the 
Luttinger Ward functional, the leading order effect of fi- 
nite temperature is taken into account by making a vari- 
ation g ^ g + Sg where 



sq = +eg, 



in other words, 



VII. LOW TEMPERATURE FREE ENERGY 

The low temperature Free energy. Recapitulation of 
Luttinger Ward's hard argument. If we use the zero 



Y[g,T\-Y[g] = 

where Y = Y[0,6\. 



dw , 
2ttI 



-Tr 



5Y 
Sy (u) 



(63) 
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It is instructive to briefly digress to see how this replacement works in the leading order contribution to the Luttinger 
Ward functional, given by 



Yi[g,T\ = KxNx (^=J T 2 E Gc(u„)g x (iv r - iw n )Q h (iv r 



(64) 



So we can write 

Y X [G,T] = KV 2 



(2m) 

To leading order in e, we have 



^2__± 2m8(iu)i + iuj 2 - iu) 3 )[l + e(u>i)][l + e(w 2 )][l + e(uj 3 )]g c (uii)g x (ui 2 )Gb(uJ3) ■ 



5Y X = KV 2 



(2m) 3 



2mS(itJi + iuj 2 - mj 3 ) [e(u>i) + e(w 2 ) + e(wa)] ^(wi)^^)^^) 



dio 

2~TTi 

dus 
2^i 



sry fry 
NK2wi-^-e(uj)g c (cj) + K2m'-e(uj)g x (u;) + N2m-^-e(cj)g b (uj) 
dy c (w) oy x (uj) oy b {uj) 



-Tr 



6g(<jj) 



(65) 



where the factors of N, NK and K arise because the functional derivative refers to the variation of Y\ with respect 
to a given spin/flavor leg of the propagator. The same series of manipulations can be carried out on diagrams of 
arbitrary order. 



Returning from the digression, let us now use the re- 
lation 



SY = - Str [ESg] 
to identify 



|^Str[E(w)<f0(w)] 



SY ~, . 

2m — = = — r>T,(uj). 

sg(uj) 

With this identification, we can rewrite (f6"3")l as 
dw 



(66) 



(67) 



Ay 



2m 



-Str 



= -T 



Strjs (5} 



Y,(uj)e(Lj)g(uj) 
T 

T=o 



(68) 



We see that the finite temperature change in the generat- 
ing functional correction term exactly cancels the change 
from the second, self-energy term in the Free energy func- 
tional, so that for small finite temperatures 



AF = 



TStrjh^-Cr 1 ) 
TStrjln^- 1 )} 



Ay 



T 
T=0 



T=0 



(69) 



as promised. 

Now in principle, this expression involves finite tem- 
perature contributions from all fields, including the frac- 
tionalized excitations. If we carry out the Matsubara 
summations, we see that 



AF 



TV J_^_ r 

3 27r dw 



Str Im In (-tr 1 



uj=-M+ 



(70) 



and the density of states of each field contributes to the 
specific heat coefficient 7 = —d 2 F/dT 2 . It is this obser- 
vation that primarily motivates our gap hypothesis, for 
with a gap in the holon/spinon spectrum, the only re- 
maining contribution to the linear specific heat derives 
from the conduction electrons 



_ 2tt 2 d 1 
3 dio 2n 



Im In (£-7(0!? 



J bj=-m+ 



note that as was previously done for susceptibilities we 
subtract the effect of the empty conduction band to get 
the impurity contribution to the specific heat coefficient. 
We have removed the tilde sign g — > g as from now on 
we restrict our derivation to the Green's function at zero 
temperature. 

In the single impurity model, the low energy thermo- 
dynamics can then be entirely expressed in terms of the 
conduction electron phase shift 

(72) 



<5 c = Imln[l-£ c °£ c ] , 



in terms of which 



7i 



E 



dS c 
du 



u=-iO+ 



The linear coefficient of the specific heat is then 



7 



7T" 

y 



E 

a v 



1 ds c 

—lmy c — - 

ir dio 



(73) 



(74) 



and we see that the renormalization of the density of 
states is entirely encoded in the frequency dependence of 
the conduction electron self energy at the Fermi surface. 
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The derivative of the self energies with respect to the 
frequencies can be calculated in a similar way to the 
derivatives with respect to the helds B A . The analogous 
expression is 



E 



div' 

—r aa , (uj, J)x] a ,Q % a , (u/)7a' (uj') 
2m 



^r aa /(w,o)?7 a ' 



lmGa> 



.(75) 



-i0+ 



The additional second term results from the discontinu- 
ous jump across the real axis: 5Q = — 2iIm£/(0 — irj). In 
the first term the contribution from the local conduction 
electrons to the summation vanishes in the band limit. 
On the other hand in the second term only the conduc- 
tion electrons contribute since the holon and boson are 
gapped and their spectral function vanishes at the Fermi 
level. Explicitly 



y^r cctifcct /(o,cj)g^(a;) 

. a' 

^iW aV '(0,0) .(76) 

— i0+ „',./ 







r du { 




-»0+ ■> 


' 2ri | 



ImGc 






The next step in our derivation is to compute the Ward 
identities that relate the frequency and field dependences 
of the electron self energies to the electron Landau pa- 
rameters. 

The derivation of the Friedel sum rules in Sec. 2 re- 
lied on the Ward identities of Eq. ([29]) which were derived 
by shifting the frequencies in diagrams of the Luttinger 
Ward functional Y[G] on by a shift 8u) A — Suq A propor- 
tional to the conserved charge Q A carried by the corre- 
sponding particle: G c (tu) — > G c (u) + Sui A ). In the Lut- 
tinger Ward diagrams such shifts correspond to shifts 
along along closed loops. Since the interaction vertex 
conserves the charges Q A it is also possible to think of 
similar shifts in the self energy diagram o 12 ' 13 ' 29 ^. The 
interesting case is when the external legs of the self en- 
ergy carry the relevant charge q A ^ 0, in this case in 
addition to the closed loops also the external frequency 
is shifted and the self energy change can be related to the 
frequency derivative of the self energy 

8uJ A ~ q A dw' [ ' 



These identities introduce in addition to Eq. (|75|) alter- 
native equations for the frequency derivative of the self 
energies 

1 dE a (u)) x-^ f du}'. 



1* 



du> 



E 



2iri 

z ^T aa ,(uj,0)q£ri a 



(u,uj')r] a >gl, (uj')~f a >q A 



(78) 



-i0+ 



The conduction electrons carry the three charges q A ^ 0, 
and it is possible to plug these derivatives into the ex- 
pression for the specific heat coefficient of Eq. ([71)1 , which 
gives 



1_ 

To 



1 



NK 



du) 



(0)W)^ a (o;)r7 a g a 7 Q 



aua 



av^ cav,a 

(0,0)g 

a 



lvaQ a 



(79) 



with 70 = NKir 2 /3 ■ ImG c /ir the specific heat coeffi- 
cient of the bare conduction band. Rewriting the latter 
equation we can relate the specific heat coefficient to the 
various susceptibilities 



To 



X 
Xo 



NK 



(80) 



where Xo = NKImQ c /ir are the bare Pauli susceptibili- 
ties. The latter expression can be simplified as a set of 
Ward identities relating the specific heat to the suscepti- 
bilities and the on-shell Fermi liquid interaction vertex 



To 



x s 


i 


Im£ c 


c 

Xo 


NK 




x F 


1 


IvaQc 


*r 


h NK 


IT 


x c 


1 


lvaQ c 


xr 


F NK 


7T 



E 



E -' r 



(0,0) 
,'(0,0) 



ca.iy.ca' , 



aa'vv' 



','(0,0). 



(81) 

In deriving these identities we assumed a gapped spinon 
and holon spectra which restrict the summation over a 
in Eq. ([80|) to the conduction electrons. As discussed in 
the following section it is possible to directly relate the 
on-shell conduction electron interaction vertex r cai , iCQ <„/ 
to the Fermi liquid amplitudes 4> av ,a'v' by comparing the 
latter Ward identities to the Nozieres Fermi liquid expres- 
sions. Note that it is possible to derive three additional 
spectral sum rules for the off-shell interaction vertices by 
equating Eq.([76)l with Eq. (|78| . 



VIII. YAMADA YOSIDA RELATIONS 

As was mentioned in the introduction the phase shift 
in a Nozieres type local Fermi liquid can be written as a 



13 



function of the energy and the occupation 

S av [u),{nk>a>,v>}] =S +8'u+ ^ (j>oiU,a'u'8nk'a'v'i (82) 

k' a' v' 

where S' is the energy derivative of phase shift and the 
impurity Landau parameter 



4>av,a'u' — (oa/, a v \<f>\aua u ) 



(83) 



is the expectation value of the interaction between quasi 
particles. The density of the quasiparticles is po + 5' /it, 
where po is the bare density of states, which implies that 
the impurity contribution to the specific heat coefficient 
is given by 7 = jq5' /irpo. For a local Fermi liquid the im- 
purity Landau parameter can be parameterized in terms 
of two numbers <f>\ and §2 through 

4>au,a'v' = 01 (1 - 5 aa >5uv') + <p2 (<W — S aa ') ■ (84) 

The response of the conserved charges Q A to their fields 
B A can be represented in terms of the phase shifts 
through 



-Tq A 

7T ' ■ 



B A • 



(85) 



The resulting susceptibilities read 

= V* &&' '<t> a u,a'v' = -<j>\ - K<j} 2 

Xo 70 ±^ v , 

F 

X 7 
Xo 7o 



(J 

^C~— = y <j>au,a'v> = {NK — 1)01 + (N — K)<p2, 

(86) 

where the Pauli susceptibilities are Xo = NKpq/tt. By 
comparing the latter equations to the set of equations 
(f8"Tj) . it is easy to read off the relations between the Lan- 
dau parameters and the on-shell interaction vertex 



■kNK 



r cQ ,, C av(0,0). (87) 



From equations (fB"6"|) it is possible to express the Fermi 
liquid parameters in terms of the reduced susceptibilities 



X°/Xo 1 fx F /Xo , X 6 /Xo 

N 



•>2 = 



NK N + K V K 

X F /Xo ~ X S /Xo 
N + K 



which produce the following relation between the suscep- 
tibilities and specific heat coefficient 



NK2-^K^—l^ + N^—l^ + ^, 

7° N + K X ° S K + N X ° f X° c 



(89) 



generalizing the Yamada Yosida relatio n 12 ' 13 i 31 to our 
infinitc-U model with K channels. Note that also the on- 
shell interaction vertex can be parameterized in a similar 
way 

Tcav,a'u' (0,0) = r 'i(l - S aa > 5 UV <) + T 2 {5 V v' —5 aa ')t (90) 

with L12 = —TrNKcfii^/ImSc- These two Fermi liquid 
parameters control the low temperature thermodynam- 
ics, and are renormalized due to intermediate spinon and 
holon states. 

An interesting observation is that the prefactors of the 
susceptibility terms in Eq. ([59|) correspond to the three 

central charges of the SU{N) K , SU{K) N and U(l) Kac- 
Moody algebras respectively 3 ^,. In terms of the relevant 
Kac-Moody currents the free fermion problem decom- 
poses to the free current Hamiltonians H^, TL F and Ti%- 
The fraction of the degrees freedom in section A, 7^/7°, 
is given by the ratio of the central charge of section A to 
the total central charge NK—. Equation (|89[) emerges 
under the assumption that for the fully screened impu- 
rity problem the fixed point Hamiltonians have a similar 
form to the bare Hamiltonians with renormalized Fermi 
velocities and no residual degrees of freedom 34 



nit 

n A 



FP 



1 



(91) 



where Xa represents the impurity effective couplings and 
I is the size of the system. Under this assumption 
Xa/Xa = 1a h\ = 3ir\ A /l, where j A is the impurity 
contribution to the specific heat coefficient in section A. 
Summing up the impurity contribution to the specific 
heat coefficient 7 = J^a^a we recover Eq. 



7 



7a Xa 



?7° X°a 



NK 



K- 



M 2 



1 y s K 2 
XS + N- 



1 Xf , Xc 



N + K X " S N + Kxf 





Xc 



IX. SHIBA -KORRINGA RELATIONS 

The Korringa-Shiba relations relate the impurity's 
nuclear-spin - lattice relaxation time T\ with the Knight 
shift K at small magnetic fields. Alternatively it relates 
also between the impurity static susceptibility and the 
low frequency impurity spin power spectrum. In the con- 
text of our paper this relation is important as a mani- 
festation of the soft modes formation in the susceptibil- 
ity power spectrum even though the spinon spectrum is 
gapped. This relation was first deduced from experimen- 
tal data by Korringa, and only later it was proved" by 
Shiba for the finite-U Anderson model summing up the 
perturbative expansion in U to arbitrary order. Below 
we show how Shiba's arguments can be extended to the 
infinite-U Anderson model and to the Kondo model by 
summing up the perturbative expansion in 1/N. 
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The expression deduced by Korringa reads 

K 2 T{T = C, (92) 

where C is a constant. The Knight shift K describes 
the mean field response observed by the nuclear moment, 
which for low Larmor frequency is proportional to the 
static impurity susceptibility \s- On the other hand the 
relaxation rate T^ 1 is proportional to the low frequency 
power spectrum of the impurity susceptibility. Explicitly 
in terms of impurity susceptibilities K and Tj -1 are given 
by 



- -T{g Ni i N fA 



2 A2 
HF 



x+-"(n-io+) 



where Ahf is the hyperfine coupling and gjy and /lijv 
are the nuclear spin g-factor and Bohr magneton respec- 
tively and £1 is the Larmor frequency. The latter ex- 
pression with Korringa's empirical rule give immediately 
the Shiba relation between the static susceptibility and 
the spin power spectrum at low frequencies where the 
susceptibilities become isotropic 



A straightforward calculation of the diagram (|94|) gives 



tco dujduj' 



tf(^(-io+)) 2 



E 



(2vri) 



g b (ui)r baca » (ui, o) 

(<W)0 6 V)- 



(94) 



In order to obtain the relation between the power spec- 
trum of Eq. ([94)) and the static susceptibility of Eq. (f57|) we 
need to resort to the symmetry properties of the vertex 
function T. In addition to the symmetry to the exchange 
of variables which is inherited from the bare interaction 
vertex A, there is an additional symmetry of the vertex 
function as a two particle vertex. In general a two parti- 
cle vertex T agt , a 'a' is a function of the incoming spins a 
and a' and of the outgoing spins a and a' . In our model 
there is no annihilation or creation of spin therefore each 
incoming spin must identify with an outgoing spin 

r - t -t = rWfi -S /-/+r (2) ^ -,s (vs) 

1 aa,a'a' — x u aa u a' a' ~ 1 u aa' u a' a- \<ju J 

In the full interaction vertex function the incoming and 
outgoing spins of each particle are already identified a = 
a and a' = a' hence we can write 



x% 



n-»0 



(93) 



Extending Shiba's arguments to the infinite-U Ander- 
son model we see that by taking the derivative of the 
dynamic susceptibility dx" /dQ and substituting ft — ► 0, 
the only contribution to the low frequency power spec- 
trum comes from the following diagram 



T caba , (w, J) = (w, J) + T[;> (w, u/)5 aa , . (96) 

The static susceptibility of Eq. (|5?|) gets a contribution 
only from the exchange part of the vertex I^ 2 ) as the 
direct part drops out of the summation ^ aa'T. 
The resulting simplified expression is 



Xs 



NK 



dw 
2^' 



6 2 Hrf 6 )(O,a;). (97) 




The line that intersects the upper conduction electron 
line denotes the non-analytic part of the derivative of the 
green function which arises from the jump across the real 
axis dg c /du>\ na , = —1iG"(—iQ + )8(iuj). Since the boson 
and holons are gapped there are no similar contributions 
from derivatives of the boson and holon green functions. 
The analytic part of the green function derivatives pro- 
duce no imaginary contribution and should not be con- 
sidered here. The other diagrams arise from derivatives 
of bare interaction vertex A. These diagrams have more 
than one electron-hole pair hence their contribution to 
the imaginary part is super linear in Q and the deriva- 
tives vanish as fl — > 0^. 



A similar simplification can be presented for the power 
spectrum of Eq. (f94|) by noticing that 



y^oT bQ 



^ a'r ca "ba' 



" 1 be ' 



~M2) 



ft "I 1 

= a 1 



cb 



(98) 
(99) 



Taking further advantage of the symmetry to variables 
exchange we can rewrite the power spectrum as 



NK(g'jf 



7 g 2 ( w )rf (o,u/ 



The right-hand side of the latter equation can be ex- 
pressed in term of the static susceptibility of Eq. ([97]) 
which gives the Shiba relation 



xS(«-«> + ) 



NK 



Xs- 



(100) 



According to this derivation the Korringa constant C of 
Eq.® is NK/[n(g N ^ N ) 2 }. 
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This result implies that the dynamical susceptibility is 
dominated at low frequencies by intermediate electron- 
hole pair states. These intermediate states of the Fermi 
liquid dress the dynamical susceptibility to produce a 
non-vanishing power spectrum at low frequencies. While 
the spin-charge decoupled excitations remain gapped, the 
are soft magnetic modes that are carried by intermediate 
quasiparticles. 



X. DISCUSSION 

This paper was motivated by an interest in the way a 
Landau Fermi liquid can emerge within the constrained 
Hilbert space of highly correlated electron systems. In 
such highly constrained systems, we can no longer ap- 
peal to adiabaticity or the use of infinite order expan- 
sions in the strength of the interaction. By combining 
the Luttinger Ward approach with a new class of large 
N expansion for the infinite U impurity Anderson model, 
we have obtained a novel perspective on the Fermi liquid 
that may have implications for our understanding of the 
corresponding Anderson and Kondo lattice models. 

In the single impurity model, we find that the Fermi 
liquid contains low lying electronic quasi particles that 
break-up into fractionalized excitations above a certain 
gap energy A fl . The assumption that this gap holds to 
all orders in our analysis leads to a full set of conserv- 
ing relationships between dynamic and thermodynamic 
variables, but it also has qualitative implications for our 
understanding of the Landau Fermi liquid. 

Indeed, our results indicate a rather intimate relation 
between the opening of a gap in the spinon and holon 
spectra, and the formation of a heavy Fermi liquid. In 
particular, the Landau scattering amplitudes that appear 
in our treatment contain, as intermediate states, the 
virtual excitation of spinon or holon pairs, for example, 




These amplitudes are proportional to the inverse of the 
gap, and are finite so long as the gap is present. Through 
this description we showed how these virtual excitations 
give rise to the renormalization of the electron-electron 
interaction, and a whole host of conserving relationships 
known to be valid for the Fermi liquid. 

But is this really a Fermi liquid? Landau's notion of a 
Fermi liquid certainly has no place for fractional quasi- 
particle excitations, even if those particles are gapped. 
On the other hand, the fractional particles that appear 
in this theory involve fields with matrix elements that 
link states of different Q: states that corresponding to 
models with different impurity spin S — Q/2. The exci- 
tation "gap" in the spinon spectrum involves excitations 
from the ground-state of the fully screened model with 



spin 5*o = K/2, to excited states of the under-screened 
model, with S = Sq + i, 

G b (t) ~ $> + i, a\bi\Sv)eA E ^>^- E ^> 

a 

^ e -i(s s [So+|]-B s [S ])t_ / 101 -j 

where E g [So] and E g [So + 5] are f ne ground-state en- 
ergies of the fully screened and under-screened models 
respectively. Thus the appearance of the gap in the 
spinon/holon spectrum is purely a consequence of the 
stability of the fully screened ground-state. 

Physical quantities, such as the spin or charge corre- 
lation functions, or the f-spectral function involve com- 
binations of spinons or holons which do not change Q 
and do not display this gap. For example, the spin rais- 
ing operator creates a spinon/antispinon pair. Such ex- 
citations are "gauge neutral" and are contained in the 
physical Hilbert space of definite Q. These pair excita- 
tions couple to the conduction sea, converting the sharp 
gap of the infinite iV limit into a pseudogap at finite N. 
The associated time-scale in this gap t c ~ ft/A g gap can 
be interpreted as a confinement time-scale, for after this 
period the spinon-antispinon pair combines with holons 
from the sea to produce an electron-hole pair 

s a + sp — > 9 e all + ep^. (102) 
Diagrammatically, 



ct \i a \i 




+ ... 



Pn Pn 



. (103) 

(This process can be loosely compared to "quark jets" in 
a collider, where quark- antiquark pairs rapidly combine 
with virtual quarks from the vacuum, to form mesons or 
baryons.) 

The processes that confine the spinons are precisely 
those that gave rise to the gapless spin spectral func- 
tion. From our derivation of the Shiba relationship, we 
know that the decay of spinon-pairs into the electron-hole 
continuum is constrained to give rise to a spectral func- 
tion with a pseudo-gap whose zero-frequency intercept is 
determined by the relation x"{ u> )/ UJ \u=o = ^^/{NK): 




where the dots will represent the DC limit ~x"(ui)/u> 
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\ 2 tt/{NK) predicted by the Shiba and the dashed line 
represents the gapped, N — > oo limit. A detailed cal- 
culation of this conserving spectrum within the leading 
Kadanoff Baym approximation^ would be an interesting 
exercise for future work. The Fourier transform of this 
power spectrum leads to the correct l/t 2 decay of the spin 
correlation function (S(t)S(0)). Thus the confinement of 
the spinons through their coupling to the electron contin- 
uum is an integral part of the emergence of Fermi liquid 
behavior in the single impurity model. The pseudo-gap 
in this function is quite unremarkable from the Fermi 
liquid perspective, where the the offset peaks in the dy- 
namical susceptibility are a well-known signature of the 
asymmetric Kondo resonance . 

In this way, the confinement of spinons and holons to 
form a Landau Fermi liquid, appears very naturally in the 
Schwinger boson approach to the single impurity, and it 
can be seen as a direct consequence of the way the gapped 
spin fluid couples to the conduction fluid. This is a very 
different kind of confinement to that considered in the 
context of insulating spin liquids, where there is a broken 
gauge symmetry and topological tunneling effects of the 
gauge field are required to understand the confinement 
mechanisrr£i£i2£. 



A. Lattice: heavy fermions 

There are some interesting implications of this kind 
of picture for the more general lattice models. Con- 
sider for example, the (K-channel) infinite U Anderson- 
Heisenberg lattice model 



H L attice = > ; C %va + H l(j) 



<i,3> 



H m (i,j). 



Here 
Hi(3) 



V_ 

N 



E 



H.C 



describes the atomic orbitals at the different sites and 
their hybridization with the conduction band, where 
ib\,.„ = y^irct e -ik-Rj crea tes a conduction electron at 
site j. The final terms H m (i,j) describes the antifer- 
romagnetic spin interaction between nearest neighbors, 
where 

n m (i,j) = -^B} j B ij . 

and = J2 a s g n ( a )b\ a b]- a creates a singlet pair of 
bosons along the bond In the lattice model, the 

conserved charges Qj — nb(j) + n x (j) replace the no- 
double occupancy constraint at each site, and as before, 
the fully screened case corresponds to Qj = K 

The manipulations that we have carried out in this 
paper can be extended to the model, but certain caveats 
apply. One of key differences in the lattice model, is that 



in the large N limit the antiferromagnetic interactions 
induce pair condensation (—(Jjj/N)(Bij) = A«) of the 
Schwinger bosons, which mark the onset of short-range 
magnetic correlations^. The important effect of the an- 
tifcrromagnetism at a mean-field level is described by 



N ^Mpl^ ,;|,)4) 



The appearance of a spinon-pair condensate is restricted 
to regions of large Jh ^ Tk- Once pair condensation oc- 
curs, the local gauge symmetry associated with conserva- 
tion of the Qi — K is broken, replaced by a global U(l) 
symmetry (where the difference of the Qi on even and 
odd sublattices is conserved). Now, individual spinon 
and holon excitations start to propagate between sites. 
The phase of the bond-variables becomes elevated into a 
compact U(l) gauge theory^. 

We now discuss the qualitative phase diagram of this 
model, using some of the insights gained from our impu- 
rity treatment. A rigorous application of our methods to 
the lattice requires that we include the U(l) gauge fields 
into the Luttinger Ward functional, a task that is beyond 
the scope of the current discussion. However, we expect 
that many qualitative aspects of the current analysis will 
survive. 

There are two variables that tune quantum fluctua- 
tions in the lattice model, (Fig. [3] ) 

• "Spin fluctuation parameter" 1/5 = 

• "Doniach parameter" t = Tk/ Jh where Tk is the 
impurity Kondo temperature. 

The first parameter tunes the strength of the quantum 
spin fluctuations and is loosely equivalent to the effect 
of lattice frustration. By tuning the value of 1/5, we 
can tune from classical antiferromagnetism, at small 1/5 
where the Schwinger boson individually condenses to 
form an antiferromagnet, to a quantum antiferromag- 
netism at large 1/5, where a spin- liquid with a gap de- 
velops. 

The second parameter tunes the strength of the Kondo 
effect. When t is small, boson pair condensation takes 
place and the spin physics is that of an insulating anti- 
ferromagnet or spin liquid with a decoupled background 
of conduction electrons. At large t > ti ~ 1, no pair con- 
densation occurs preserving the locality of the spinon and 
holon propagators. The direct extension of our methods 
to this case predicts a "local" heavy electron state which 
develops a large Fermi surface volume and enhanced mo- 
mentum independent scattering amplitudes. In the sin- 
gle impurity model, the charge of the localized state is 
given by the difference of the electron and the holon phase 
shifts Eq.® 



K 



NK- 



,5r 
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(a) 

2-impurity 




FIG. 3: Showing the large N phase diagram of two impu- 
rity and lattice models as a function of the ratio t = Tk/Jh- 
Boson pair condensation (66) 7^ develops for t < ti in- 
ducing short-range magnetic correlations that compete with 
the "Kondo" screening, (a) Computed large N phase dia- 
gram of the two impurity Kondo problem^ displaying the 
relative temperature scales. T — Tk and T = Tf denotes the 
lower and upper scales associated with the two stage screen- 
ing process^. The Fermi temperature and spinon gap close at 
the Varma Jones quantum critical point at t = ti, which sep- 
arates the valence bond state from the two impurity Kondo 
state, (b) Speculative sketch of the zero temperature phase 
diagram for the Kondo lattice, based on the discussion in the 
text, showing dependence on the spin fluctuation parameter 
1/5 = N/(2S). At large 1/5 and small t, a spin liquid (SL) 
forms. At small 1/5 and small t — Tk/Jh, the Schwinger 
boson pair condenses to produce an antiferromagnet. The 
spin-liquid/ AFM and spin-liquid/HF^^ phase transition are 
expected to merge into a single quantum critical point that 
separates the heavy electron state from the antiferromagnet. 
The spinon and holon gaps close at each of these quantum 
phase transitions. 

The corresponding Luttinger theorem for the latticed 
tells us that the charge density is given by the differ- 
ence of the conduction electron and holon Fermi surface 
volumes 



This expression holds some interest for us, because in 
the spin liquid phases of the model, where the Kondo 



effect is inactive, we expect the Fermi surface volume of 
the holons to be zero, giving a small Fermi surface for 
the electrons. By contrast, at large t = Tk/Jh, where 
the Kondo effect occurs, the Fermi surface volume of the 
gapped holons is now at a maximum v x / (2tt) d = 1, caus- 
ing the electron Fermi surface to expand by an amount 
Avps = l/N, forming a large Fermi surface. So long as 
the holons preserve their gap, the volume of these Fermi 
surfaces will remain fixed. 

The effect of tuning the ratio t = Tk/Jh has already 
been considered in the context of a two impurity Kondo 
model 10 (Fig. 3. (a)). Fig. 3 (b) shows the schematic 
behavior expected in the large N lattice model. Here, 
depending on the size of 1/S, either a spin liquid (large 
1/5) or an antiferromagnet (small 1/S) develops at small 
Kondo coupling t = Tk/Jh- At large t however, the 
heavy fermi liquid, with a gap to spinons and holons will 
remain stable, with Fermi liquid scattering amplitudes 
that are determined by the exchange of virtual spinons 
and holons. Provided the phase transition between these 
two limits is second order, then the spinon-holon gap 
must close at the transition line, with the concomitant 
divergence of the Landau scattering amplitudes. 

The closing of the spinon gap may have important 
implications for the evolution of the Fermi liquid pa- 
rameters near the antiferromagnetic instability. In the 
Schwinger boson scheme, at the antiferromagnetic in- 
stability, the spinons condense at one half the magnetic 
wave-vector, becoming gapless at two points in momen- 
tum space q = ±Qo/2. This will give rise to strong elec- 
tron -electron scattering at 



<f=Qo/2±Qo/2 = 




(105) 



as illustrated in 

q~Q/2 




q~±Q/2 

The appearance of strong scattering at a staggered 
wave-vector is well-known in a spin-density wave sce- 
nario. However, the additional appearance of strong scat- 
tering at q = 0, is a new feature. In other words, a sub- 
leading effect of the antiferromagnetic instability driven 
by spinon condensation, is the development of strong for- 
ward scattering reminiscent of an almost ferromagnetic 
metal. Features of this sort have been observed near the 
field-induced quantum critical point in YbRh,2Si2 39 ■ 

These qualitative features of the large N limit of our 
model are a topic of current active examination, and we 
hope to report on them in greater detail in forthcoming 
work. 
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